Microscopic dynamics underlying the anomalous diffusion 
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The time dependent Tsallis statistical distribution describ- 
ing anomalous diffusion is usually obtained in the literature 
as the solution of a non-linear Fokker-Planck (FP) equation 
[A.R. Plastino and A. Plastino, Physica A, 222, 347 (1995)]. 
The scope of the present paper is twofold. Firstly we show 
that this distribution can be obtained also as solution of the 
non-linear porous media equation. Secondly we prove that the 
time dependent Tsallis distribution can be obtained also as so- 
lution of a linear FP equation [G. Kaniadakis and P. Quarati, 
Physica A, 237, 229 (1997)] with coefficients depending on 
the velocity, that describes a generalized Brownian motion. 
This linear FP equation is shown to arise from a microscopic 
dynamics governed by a standard Langevin equation in pres- 
ence of multiplicative noise. 

PACS number(s): 05.10.Gg , 05.20.-y 



I. INTRODUCTION 

Recently, the Tsallis thermostatistics has received con- 
siderable attention in light of its growing application to a 
variety of physical systems [1] . The research has focused 
both on fundamental and phenomenological aspect of the 
issue [2]. 

Particular attention has been devoted to the issue of 
anomalous diffusion, where a significant experimental ev- 
idence has been gathered (see Ref. [3] for a detailed bib- 
liography). The description of a diffusive process (either 
classic or anomalous) is performed generally by adopting 
a time dependent formalism. The Tsallis distribution 
namely 
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using the maximum entropy principle under the con- 
straint of conservation of particle number and energy, by 
solving the variational problem: 6 [Sq — PE — aN] = 0. 

Similarly to the classic Boltzmann distribution, the 
Tsallis distribution can be also obtained as the steady- 
state distribution of a time dependent Fokker-Planck 



(FP) equation. Recently, the research on the derivation 
of the Tsallis distribution from FP equations has pro- 
duced considerable results [3-13]. The research in this 
area can be classified in one of two classes. 

First, linear FP equations are considered with diffu- 
sion and drift coefficients depending on the velocity. The 
dependence is chosen to lead to the Tsallis distribution 
as the equilibrium solution of the FP equation. Within 
the linear approach, two different choices of the drift 
and diffusion coefficients have been proposed. Stariolo 
[4] chooses a constant diffusion coefficients and alters the 
drift coefficient to include a generalized potential depend- 
ing on the Tsallis parameter q. This approach introduces 
a more relevant modification of the classic Browninan ap- 
proach. In ref. [5], instead, the classic Brownian drift 
coefficient have been considered, but with a modified dif- 
fusion coefficient to include a quadratic velocity depen- 
dence. The two linear approaches described above are 
in reality just examples of an infinite class of linear FP 
models that give a Tsallis equilibrium distribution [6]. 
Clearly, the selection of a specific linear model among 
the class requires the introduction of other criteria be- 
yond the simple requirement of leading to an equilibrium 
Tsallis distribution. 

Second, non-linear FP have been shown to lead to equi- 
librium Tsallis distributions. This approach, introduced 
by Plastino and Plastino [7] and continued by various au- 
thors [3,8-13], introduces a diffusion coefficient depending 
on powers of the distribution function. The drift, instead, 
can be equal to zero or described as in the classic Brow- 
nian motion. This latter approach, besides its elegance 
and simplicity, admits time dependent solutions charac- 
terized by retaining at every time the form of a Tsallis 
distribution. This self-similarity of the evolution repre- 
sent an important property of the non-linear approach. 

The present paper deals with the question of whether 
the linear [5] and non-linear [7] FP approach to the 
derivation of the Tsallis distributions are equivalent. The 
answer proved here is that indeed the two approaches are 
equivalent, in the sense that they both allow the presence 
of self-similar transients where the system is character- 
ized by the Tsallis distribution at every instant. 

In order to explain the microscopic origin of the 
anomalous diffusion associated to the non-linear FP 
equation of ref. [7] Borland suggested a feedback from 
the macroscopic level to the microscopic one [13]. In the 
present work, we show that the non-linear FP equation 
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of rcf. [7] but also the well known in the literature non- 
linear porous media equation, considered recently in the 
frame of Tsallis thermostatistics in ref. [8], can be recast 
in the equivalent linear FP equation of rcf. [5] . This im- 
portant result allows a deeper interpretation of the non- 
linear FP equation describing the anomalous diffusion in 
terms of a linear Langevin microdynamics in presence of 
a multiplicative noise. 

The present work is organized as follows. In Section II, 
a Generalized Brownian (GB) motion is derived from the 
Langevin equation in presence of multiplicative noise. In 
Section III, the GB motion is shown to lead to a macro- 
scopic motion described by the linear FP equation of ref. 
[5] that admits as a solution a class of time dependent 
Tsallis statistical distributions. In Section IV and V, the 
same distribution are shown to represents states governed 
also by the non linear FP equation of ref. [7] and by the 
non-linear porous media equation respectively. Finally, 
in Section VI conclusions are drawn. 



II. GENERALIZED BROWNIAN MOTION 

We consider the microscopic process described by the 
following Langevin equation: 
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The quantity —mh{t, v) is the deterministic force act- 
ing on a particle of mass m and velocity v{t) while 
mg{t,v)r{t) is a stochastic force acting on the parti- 
cle, with T{t) a Gaussian random variable with zero 
mean and delta correlation function. The presence of 
g{t,v) in eq.(l) implies that the particle is subject to a 
multiplicative noise. The distinction between additive 
(when g{t,v) = const.) and multiplicative noise (when 
g{t, v) ^ const.) is very significant when g{t, v) is a time 
dependent function. In this case arise naturally the ques- 
tion related to the definition of the stochastic integral 
(Ito or Stratonovich definition). For a more detailed dis- 
cussion on multiplicative noise, sec rcf. [14]. The micro- 
scopic process described by eq. (1) implies a macroscopic 
process described by the following linear FP equation: 
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where the drift coefficient J{t, v) and the diffusion coef- 
ficient D[t, v) have the following expression 



J{t,v) = h{t,v) 



D{t,v)=g{t,v)^ 



(5) 



(6) 



obtained using the Ito definition for the stochastic inte- 
gral. Note that, for Brownian motion: 



J{t,v)=j{t)v , 

D{t,v)=c{t) , 
the drift current in eq.(4) 
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is simplified as: 
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and the current velocity jdri ft /p becomes simply propor- 
tional to the viscous force —mhit, v) = —mjit)v of the 
microscopic process. 

A problem arises in conjunction to the results just ob- 
tained. Whether other motions, besides the Brownian 
motion, are characterized by a current velocity propor- 
tional to the viscous force. This issue corresponds to the 
existence of other solutions of the following equation for 
the unknown functions £>(t, v) and J(t, v) 

At,v) + ^^^^ = eit)Jit,v) , (11) 

in addition to the solution (7), (8), relative to the Brow- 
nian motion. The issue is easily resolved and other so- 
lutions can be founded. The more general solution is 
formed by copies of functions J(f, v) and Dit, v) where 
Dit,v) is given by 

Dit, v) = cit) + [eit) Jiv)dv , (12) 

while Jit,v) remains arbitrary. The simplest solution, 
for which J(f, v) is given by (7) provides the definition 
for a new generalized Brownian (GB) motion [5]. 



III. LINEAR FOKKER-PLANCK EQUATION 

We consider the FP equation (4) for the GB processes. 
With the introduction of the dimensionless time r: 



dr = eit)^it)dt , 
and the functions -D(r), /3(r) and parameter q: 
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Prom eq. (19) it follows 

Z,(r)/3(r)V2 = Z,(0)/3(0)V2 = TV, . (24) 

The constant Nq is determined starting from the expres- 
sion of Zq{T) given by 



the diffusion coefficient (12) with drift coefficient given 
by (7) can be written in the following form 
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D{t, v) = D{t) [1 - (1 - q)(i{T)v^] , (16) For g > 1 [15], it results: 

while after taking into account (7), the FP equation (4) 
becomes [5]: 
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The final solution of eq. (17) has the form 

p(r, v) = iV,/3(r) [l _ (l _ q)piry] ''^'-'^^ 
where (3{t) is given by eqs. (22-23). 
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The time-dependent solutions of eq.(17) are sought us- 
ing the following ansatz: 
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The above ansatz requires the solution to conserve at 
every time the form of a Tsallis distribution with time 
dependent parameters Zq and /3. The time dependence 
of the two parameters determines the actual solution and 
is obtained easily substituting ansatz (18) in eq.(17). It 
follows that the equations determining the evolution of 
Zq{T) and P{t) are identical to the equations for the 
Brownian motion: 
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The result above justifies the use of the term Generalized 
Brownian motion used to name the process defined by 
eqs. (7,12). Prom eq. (20) the condition below follows 



2/3(oo)£)(oo) = 1 
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again in complete similarity with Browinan motion. Eq. 
(20) is solved easily with the the substitution y = (3~^ 
that linearizes the equation: 
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IV. NON LINEAR FOKKER-PLANCK 
EQUATION 

The scope of the present and the next section is to show 
that the time dependent solution (27) obtained here, of 
the linear FP equation (17) proposed in [5], is also solu- 
tion of non linear FP equations which can be obtained 
from the linear FP (17). The goal of the present section 
is to investigate the relationship between the linear FP 
(17) and the non-linear FP equation proposed by Plas- 
tino and Plastino [7]. 

We start the proof by noting that eq. (27) allows us 
to write 

1 - (1 - qmry = 7V«-i/3(r)(«-i)/2p(T,t;)i-« . (28) 
Besides the following time dependent function is defined: 
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Then, it follows that eq.(17) can be rewritten as: 
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that is identical to the equation proposed by Plastino and 
Plastino that was solved using the same ansatz (18) used 
above but assuming that -Di(t) is constant. 

The procedure outlined in the previous section leads to 
the same relationship (19) between Zq{T) and /3(t) and 
/3(t) is governed by the following evolution equation: 

^ = 2/3(t) - 2/3(oo)(«-^)/2^1^lI/3(r)(5-«)/2. (31) 
dr Di[(X)) 

The condition (21) transforms now as 
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2/3(oo)(3-«)/2i?,(oo) = -^ 
2-q 



(32) 



As above, a transformation y = /j'^ ^'/^ linearizes 
eq.(31) and the general solution follows easily: 



/3(t) = /3(oo) 1 + 
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that is solved readily to obtain 
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From eq. (37) and (39) it follows that Z?2(oo) = 0, a 
condition that guarantees that the current (38) vanishes 
for T — > oo, as needed to obtain an equilibrium state. 
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The complete solution of eq. (30) is given by eq. (27) 
where now /3(t) is expressed as a function of £>i(t) by 
eqs. (33-34). In the special case of Di{t) constant, the 
results presented in the literature [3,7] are recovered. 



V. NON LINEAR POROUS MEDIA EQUATION 

For the solutions (27) considered above, the following 
results follows readily: 
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Substituting the relationship above in eq.(30) it follows: 

(36) 



where 



D2{t) - ^/3(r)(«-3)/2 
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As a consequence, eq. (36) is the well known non-linear 
porous media equation, widely used in condensed matter 
physics and considered recently in ref. [8]. The current 
of particles is given by: 



j-(r,z;) = -i?2(r)|-b(r,t;)]2-« 



which generalizes the Fick law (indeed, for q = 1 the 
classic Fick law is recovered). 

Clearly, the time dependent Tsallis distribution (18) is 
solution of eq. (36) as well. Zqir) and /3(t) are connected 
by eq. (19) and the solution acquires the form (27) while 
the evolution law for /3(t) follows from D2{t). The final 
differential equation is 



^=4(g-2)7Vi-«i^,(r)^(r)(5-e)/2 , 



VI. CONCLUSIONS 

In the present work, the same time dependent Tsallis 
statistical distribution given in eq. (27) is solution of all 
three equations considered, eq. (17), eq.(30) and eq. (36) 
and describes anomalous diffusion. The time evolution of 
(3{t) is the same in all the cases considered and can be 
expressed in terms of the three functions D{t), Di(r), 
D2{t) that are inter-related. 

The relationships derived above allow us to interpret 
the states (27) at a microscopic level. Such states are so- 
lutions of non linear FP equations (eqs. (30) and (36)) and 
consequently describe anomalous diffusion (when q 1). 
But the results obtained above, prove that anomalous 
diffusion can also be described with the linear FP equa- 
tion (17), which has a variable diffusion coefficient given 
by eq. (16). 

The primary results of this equivalence is that the 
linear FP equation (17) can be related directly to the 
microscopic dynamic model expressed by the Langevin 
equation (1) linking in this way macroscopic processes 
described by anomalous diffusion with microscopic pro- 
cesses characterized by multiplicative noise. 
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